Kink-Antikink Formation from an Oscillation Mode by Sudden Distortion of the 

Evolution Potential 
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We demonstrate numerically that an oscillation mode in 1+1 dimensions (eg a breather or an 
oscillon) can decay into a kink-antikink pair by a sudden distortion of the evolution potential which 
occurs within a certain time or space domain. In particular, we consider the transition of a sine- 
Gordon potential into a $ 4 potential. The breather field configuration is assumed to initially evolve 
in a sine-Gordon potential with velocity v and oscillation frequency uj. We then consider two types 
of numerical experiments: a. An abrupt transition of the potential to a $ 4 form at to = over 
the whole 1-dimensional lattice and b. The impact of the breather on a region x > xo = where 
the potential has the $ 4 form which is different from the sine-Gordon form valid at x < xo = 0. 
We find that in both cases there is a region of parameters (v,ui) such that the breather decays to 
a kink-antikink pair. This region of parameters for kink-antikink formation is qualitatively similar 
with the parameter region where the energy of the breather exceeds the energy of the kink-antikink 
pair in the <3? 4 potential. We demonstrate that the same mechanism for soliton formation is realized 
when using a gaussian oscillator (oscillon) instead of a breather. We briefly discuss the implications 
of our results for realistic experiments as well as their extension to soliton formation in two and 
three space dimensions. 

PACS numbers: 



1. INTRODUCTION 

Topological solitons (defects) [l|, 0, 9 are static lo- 
calized solutions of nonlinear partial differential equa- 
tions which are stable due to nontrivial topological prop- 
erties of the vacuum manifold. These coherent non- 
perturbative excitations are distinct from the perturba- 
tive 'particle' excitations which correspond to localized 
small field oscillations around the vacuum [4(. 

Solitons arise in the context of effective field theories 
[H> @ (eg Ginzburg-Landau theory) and play an impor- 
tant role in many branches of modern physics. For ex- 
ample they form in condensed matter systems Q (liq- 
uid Helium, liquid crystals or superconductors) and in 
cosmology @, E3, EH (cosmic strings, domain walls, 
monopoles) in the context of high energy physics (HEP) 
models. Solitons usually form during phase transitions in 
condensed matter or cosmological systems. In the con- 
text of HEP models, these phase transitions are only re- 
alized in the early universe and it is currently not possible 
to reproduce them in the laboratory or in accelerators. 
Thus, the experimental study of the formation of topolog- 
ical defects predicted in HEP models can only be made by 
drawing analogies with condensed matter systems where 
the experimental realization of such phase transitions is 
possible 0,0 • 

The understanding of the transition between the 'par- 
ticle' and the 'soliton' sectors is an important unsolved 
problem. The development of methods that would allow 
the formation of solitons from particles either through 
scattering or through decay would open new possibilities 
for the experimental study of solitons predicted by HEP 



theories. 

Previous approaches to this problem have focused on 
the formation of solitons from scattering of two [HI, [l3| 
or more [4] particles. These approaches are motivated by 
a time reversal of a soliton-antisoliton annihilation which 
is accompanied by emission of oscillation modes (parti- 
cles). In the case of two particle scattering, it has been 
shown that the soliton formation is exponentially sup- 
pressed [l2| even though it may proceed more efficiently 
in the presence of a pre-existing soliton 0, [lj| . A multi- 
particle scattering has been shown to be more effective in 
producing solitons [H through resonant build-up of par- 
ticle oscillations which allows the oscillations to extend 
to other vacua. The initial conditions for this process 
however require careful tuning in order to achieve the 
required resonance [3| • 

An alternative approach to the problem of soliton for- 
mation involves the decay of a single highly energetic os- 
cillation mode (particle) to a soliton-antisoliton pair (eg 
kink-antikink in 1+1 dimensions). This decay may be 
facilitated by a time or space dependent distortion of the 
potential determining the dynamical evolution of the os- 
cillation mode. Such a distortion may be achieved either 
by changing an external parameter of a system (pressure, 
temperature, external field etc) or, in the space depen- 
dent case, by considering two different materials sepa- 
rated by a given surface. 

The type of distortion considered in the present study 
is distinct from the corresponding distortion occurring in 
the effective potential during a phase transition. First, 
our initial conditions are not thermal (they consist of 
an oscillator mode) and second, we only consider mild 
forms of distortions that do not change the topological 
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FIG. 1: a: The energy of the breather in the sine-Gordon Lagrangian is smaller than the energy of a kink-antikink pair for all 
values of uo (see also eq. (|2.7p ). Here we used v = 0.3 but a similar plot is obtained for all values of v. b: The relative energy 
difference of the breather in the sine-Gordon Lagrangian with respect to the sine-Gordon kink-antikink pair is independent of 
the boost velocity since both field configurations (breather and kinks) are exact solutions of the sine-Gordon Lagrangian. 



properties of the vacuum. In particular, we consider the 
transition from a sine-Gordon potential to a <E> 4 potential 
in 1 + 1 dimensions occuring at a given time t — to or a 
given spatial point x — xq. The initial field configuration 
is an oscillation mode in the form of either a sine-Gordon 
breather or a gaussian oscillator ('oscillon'). The energy 
of the initial field configuration as recorded by the dis- 
torted potential is tuned by the boost velocity, the field 
configuration parameters (eg the oscillation frequency) 
and the phase of the configuration at the space-time point 
of the transition (ie the value of t or of xq). 

An important condition required for the decay of the 
initial oscillation mode to a kink-antikink pair is the fact 
that the energy of the oscillation mode as recorded by 
the distorted potential ($ ) should be larger than the 
energy of the kink-antikink pair in the same potential. 
This necessary conditions is of the form 



• A gaussian oscillator (oscillon) 

We compare the theoretically allowed (based on energet- 
ics) parameter region of decay with the corresponding 
region based on numerical simulations of field evolution 
and we find qualitative agreement. We also demonstrate 
numerically that the decay of the boosted sine-Gordon 
breather in the context of the sine-Gordon Lagrangian 
(without switch to <& 4 ) does not occur for any boost ve- 
locity as expected by energetic arguments and Lotentz 
invariance since 
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where E® m is the energy of the boosted oscillation mode 
evaluated in the distorted potential (<I> 4 ) at the transition 

time to and E kink (v = 0) is the corresponding energy of 
the kinks in the same potential. We have considered 
the conservative assumption that the kink-antikink pair 
forms at velocity. In practice, momentum conservation 
implies that the kink-antikink velocity is non zero. This 
is expected to decrease somewhat the allowed parameter 
region for soliton formation. 

The goal of the present paper is to investigate the decay 
of a particle oscillation mode into kink-antikink pairs in 
the context of a <E> 4 potential in 1 + 1 dimensions. We 
investigate the v—co (boost velocity-oscillation frequency) 
parameter region where the decay occurs and consider 
two types of particle-like oscillators 

• A sine-Gordon breather 



2. TIME DEPENDENT DISTORTIONS OF THE 
POTENTIAL 

As discussed in the Introduction, particles at the clas- 
sical level may be represented as small field oscillations 
around the vacuum. At the quantum level, these oscilla- 
tions are quantized and the energy of the lowest quantum 
state is equal to that of a particle. Thus, the lowest en- 
ergy oscillation mode may be identified with the particle 
excitation of the theory. A useful field configuration that 
can play the role of an oscillation mode at the classical 
level is the breather solution of the sine-Gordon model[l|. 
The sine-Gordon Lagrangian is of the form 



i 7T 



with field equation 
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FIG. 2: The boosted breather evolved by the sine-Gordon 
Lagrangian for (v,u>) = (0.5,0.6). As expected, the decay to 
a kink-antikink pair does not occur for this and for any other 
parameter value we have tested. 



The breather solution is given by [l[ 
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<& e b °(t, x; u,v) = -1 + - tan- 1 
where 

T = j[t — v(x — xq)] , X — j[x — xq — vt] 

7 = (i-« 2 r 1/2 , v 
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In addition to the breather, the sine-Gordon Lagrangian 
has soliton solutions (kinks) that interpolate between 
neighboring vacua = 2n + 1, $ n +i = 2n + 3, (n = 
—1,0, +1, ...). Such a sine-Gordon kink interpolating 
between the vacua $_i = —1 and $q = +1 is of the form 



0jEg, fc (t,a;;w,w) = -l+-tan 1 [exp(X)] 
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is straightforward to show that the sine-Gordon energy 
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For the breather Q2.3|) , the energy (|2.6|) is smaller than 
the corresponding energy of the kink-antikink pair (|2.5D 
ie 
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for all values of w, u>. This is shown in Fig. 1. We have 
verified numerically that ^(v,w) is independent of the 
boost velocity since both the breather and the kink are 
exact solutions of the sine-Gordon field equation (see Fig. 



lb) and therefore E§ g (v,oj) 



1 E b {0,u) aadEf£ k (v) = 
7-Efcmfc(0). We have verified the stability of the sine- 
Gordon breather in the context of the Lagrangian (|2.1|) 
by numerically simulating its evolution for various values 
of v, to. Snapshots of such a simulation are shown in Fig. 
2. 

We now consider the embedding of the breather as an 
initial condition in a $ Lagrangian of the form 
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We will consider a numerical experiment such that the 
breather evolves in a sine-Gordon Lagrangian (|2.ip where 
it is an exact solution for t < but at t = t = the 
Lagrangian changes abruptly to (|2.8|) . Since the evolu- 
tion for t < is trivial, we only consider the evolution 
for t > 0. The equation of motion for t > is 

$ = $" - ($ 2 - 1)$ (2.9) 

and the corresponding energy and momentum are 
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In addition to perturbative particle excitations around 
the vacua this equation of motion admits static non- 
perturbative kink solutions of the form 



P=- dx§ 
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which interpolate between the two vacua <f> = ±1. One 
way to mimic the presence of particles in the $ 4 poten- 
tial is to use arbitrary field configurations that oscillate 
around a single vacuum with long lifetimes. One such 
choice is the sine-Gordon breather (|2.3p embedded in the 
$ Lagrangian. 

The sign of ^M- depends not only on»,w but also on 
the phase of the breather at the time to of the potential 



switch. Thus 
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{v,u),t(s). In what follows we fix to 



(to = 0) and we focus on the parameters v, uj demonstrat- 
ing that kink-antikink formation is possible for a range 
of these parameters. 

The range of boost velocities v and corresponding val- 
ues of uj where energetics allow the decay of the boosted 
breather to a kink-antikink pair is shown in Fig. 3 
(shaded region in v — u> plane) . Even though, energetics 
allow the decay of the embedded sine-Gordon breather 
into a kink-antikink pair for a range of parameters, it is 
not necessary that this decay will occur in practice. For 
example the extra energy may be emitted in the form 
of radiation or remain as breather oscillating energy. In 
order to find if the decay actually occurs we have per- 
formed numerical simulations of the embedded breather 
evolution spanning the full range of the v — u> parameters. 
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FIG. 3: The parameter range where ^ > is indicated by 
the shaded region. The dots indicate parameter values where 
the decay of the boosted breather into a kink-antikink pair 
was observed in the simulations. The whole parameter region 
shown in the plot was scanned at steps of Av — 0.05, Alj — 
0.05. 



We have performed the simulations using the ND- 
Solve routine of Mathematics [2(| . solving the partial 
differential equation (|2.9[) with the embedded breather 
(|2.3p as initial condition and fixed boundary conditions 
<J> = — 1, $ = 0. Using these simulations we have con- 
firmed that there is a parameter range where the embed- 
ded breather spontaneously decays into a kink-antikink 
pair at the time of the potential switch. In Fig. 4 
we show the evolution of the embedded breather field 
for two parameter sets v — lu corresponding to stability 
(Fig. &a,:(v,to) = (0.5,0.8)) and instability towards kink- 
antikink decay (Fig. 4b:(u,w) = (0.7,0.3)). 

We have performed various tests to secure the validity 
of our numerical simulations. These tests include the 
following: 

• We have verified that total energy (eq. (|2.10p ) was 
conserved at a level of about 1% during the evolu- 
tion. 

• We have verified that total momentum (eq. (|2.1ip ) 
was conserved at a level of about 1% during the 
evolution. 

• We have doubled the size of the lattice and verified 
that the evolved field configurations remain practi- 
cally unchanged. 



• We have considered 5 different numerical algo- 
rithms for the solution of the partial differential 
equation (|2.9[) and they all lead to practically iden- 
tical field evolution. The algorithms were used 
through the 'Method' option of the NDSolve rou- 
tine of Mathematica and include the methods: 
Adams, ExplicitRungeKutta, ImplicitRungeKutta, 
DoubleStep with submethod: ExplicitModified- 
Midpoint and Extrapolation with submethod Ex- 
plicitEuler . 

The instability parameter range obtained using simula- 
tions is spanned with dots in Fig. 3 and is smaller than 
the shaded range anticipated from energetic considera- 
tions. Indeed, as discussed above, the supply of sufficient 
energy is a necessary but not a sufficient condition for 
the kink-antikink formation. For example, center of mass 
momentum needs to also be conserved for kink-antikink 
formation while in the energy criterion ^S- was de- 
fined assuming zero velocity of the pair at formation. The 
assumption of a nonzero velocity of the kink-antikink at 
formation would tend to decrease the theoretically al- 
lowed parameter region from energetics bringing it to 
even better agreement with simulations. 

The embedded oscillating breather ansatz has the ad- 
vantage of being an exact solution to a Lagrangian simi- 
lar to that of the <I> 4 potential which implies a long life- 
time and minimal radiation emission during its evolu- 
tion. On the other hand, it has a maximum angular 
frequency Lu max — 1 and this limits the maximum energy 
that can be achieved by boosting. As a result we have 
not observed more complicated decay products than a 
kink-antikink pair (eg two kink-antikink pairs) since even 
though we can easily obtain > 2Mki n k it is not pos- 
sible to achieve E®. L > 4Alkink which would allow decay 
to two kink-antikink pairs. 

In order to bypass this limitation we have considered 
an alternative oscillating ansatz corresponding to a gaus- 
sian oscillator also known as 'oscillon' in previous studies 
[3 HH [H [H, H]. This gaussian oscillator field config- 
uration is of the form 

& go {t, x; ,v,u>) — -l + asm(cuT)exp(-bX 2 ) (2.13) 

where a, b are parameters that we fix to the values a = 1 
and b = 0.1 in most of what follows. 

The v — lu parameter region where > and there- 
fore the decay is energetically allowed, is shown in Fig. 
5 (shaded region). The region on the right of the dashed 
line corresponds to parameter values where the gaussian 
oscillator energy is larger than two kink-antikink pairs 
(four solitons) implying the possibility of more compli- 
cated decay products. 

In order to confirm the above expectations for boosted 
oscillator decay we have performed numerical simulations 
for the evolution of the boosted oscillator of equation 
(|2.I3[) in a <I> 4 potential for various v — lu parameter val- 
ues. We have confirmed the existence of parameter re- 
gions where the decay to a kink-antikink pair occurs and 



5 



v=0.5, w=0.8: No Breather Decay 
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FIG. 4: a: The field evolution of the boosted breather with (v,ui) = (0.5,0.8), embedded in the <1> 4 Lagrangian. The breather 
evolves oscillating without decay to solitons since there is not enough energy for the decay, b: The field evolution of the boosted 
breather with (v,ui) = (0.7,0.3), embedded in the $ 4 La grangian. In this case, the boost provides sufficient energy for the 
decay to a kink-antikink pair. 
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FIG. 5: The parameter range where W- > is indicated by the 
shaded region. The dots indicate parameter values where the 
decay of the gaussian oscillator into a kink-antikink pair was 
observed in the simulations. Thicker dots indicate decay into 
more complicated products (eg more than one kink-antikink 
pair or additional particles) . The energy of the initial gaussian 
oscillator is larger than two kink-antikink pairs on the right 
of the dashed line. 

we have found that this region of parameters is in good 
agreement with the expectations based on energetic ar- 
guments. In Fig. 6 we show two types of boosted oscilla- 
tor evolution corresponding to stability (v, u) = (0.7, 0.5) 
and to decay (v,oj) = (0.5, 0.9). The full range of param- 



eters where we observed decay of the oscillator to a kink- 
antikink pair is spanned by small dots in Fig. 5 while 
thick dots correspond to the observation of more compli- 
cated decay products involving either two kink-antikink 
pairs or a kink-antikink pair plus particle like oscillators 
in either of the two vacua. Clearly, the region spanned 
by the dots where decay occurred in the simulations is 
in good agreement with the shaded region predicting the 
decay on the basis of energetic arguments. 

3. SPACE DEPENDENT DISTORTIONS OF 
THE POTENTIAL 

We now consider numerical experiments in 1+1 dimen- 
sions where the dynamical evolution potential is of the 
following form 

V($) = -4[1 + cos(tt$)] x<0 (3.1) 

V($) = - ($ 2 - l) 2 x > (3.2) 

ie sine- Gordon for x < and <& 4 for x > 0. Thus, the 
potential switch occurs now in space rather than in time, 
at the point xo — 0. Our initial condition consists of an 
oscillating mode (sine-Gordon breather eq. (|2.3[> ) with 
frequency w centered at x = — 50 at t — boosted with 
velocity v. The boundary conditions are fixed as in the 
previous section. The results of our numerical experi- 
ments are shown in Fig. 7 where the parameter space 
v,u> has been spanned at steps of 0.05. These results 
indicate the following: 

• There is a range of v,u> parameters (empty regions 
in Fig. 7) where the oscillator mode is not con- 
verted to a kink-antikink pair. Instead, the mode 
either continues oscillating on the <f> 4 side (x > 0) 
or breaks in two modes evolving on both sides of 
the transition point xq = 0. 
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v=0.7, tt)=0.5: Gaussian Oscillator Docs Not Decay 



-0.5, w-0.9: Gaussian Oscillator Decays to Kinks 



FIG. 6: a: The field evolution of the boosted gaussian oscillator with {v,lo) = (0.7,0.5), embedded in the $ 4 Lagrangian. The 
'particle' evolves oscillating without decay to solitons since there is not enough energy for the decay, b: The field evolution 
of the boosted gaussian oscillator with (v,u) = (0.5,0.9), embedded in the <fr 4 Lagrangian. In this case, the boost provides 
sufficient energy for the decay to a kink-antikink pair. 
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FIG. 7: The dots indicate parameter values where the decay 
of the breather into a kink-antikink pair was observed in the 
simulations at x$ = where the potential switches to the "I? 4 
form. Dots indicate kink-antikink formation from the decay of 
breather initial conditions, with the different colors indicating 
the different potentials at which the decay takes place. Blue 
dots indicate decay in the region of sine-Gordon potential, 
red dots indicate decay in the region of the $ 4 potential, and 
purple dots indicate decay in both regions. 

• There is a range of v,u> parameters where the os- 
cillator mode is converted to a kink-antikink pair 
with the two solitons evolving on both sides of the 
transition point xo = (purple dots in Fig. 7). Al- 
ternatively, the two solitons may be reflected back 
in the sine-Gordon side (x < 0, blue dots in Fig. 7) 
or evolve in the <I> 4 side (x > 0, red dots in Fig. 7). 

Representative frames from the simulations correspond- 
ing to the above types of evolutions are shown in Fig. 8 
Similar results are obtained when using a gaussian oscil- 
lator mode as an initial condition. These are shown in 



Fig. 9 



4. CONCLUSION-DISCUSSION 

We have shown using numerical simulations that 
abrupt distortions of the evolution potential in time or in 
space can lead to a conversion of an oscillation mode to 
a kink-antikink pair in 1+1 dimensions. The parameter 
range where this conversion occurs can be approximately 
obtained in the case of time-dependent distortions by us- 
ing energetic considerations. We have made no attempt 
to use such considerations in the case of space dependent 
distortions due to the complexity of the various possible 
outcomes of the numerical experiments involving multi- 
ple soliton evolution in more complex potential forms. 

Since, the oscillating modes are used to represent par- 
ticles at the classical level, the above described phe- 
nomenon may correspond to a new mechanism for soli- 
ton formation from particle decay. It may be effective 
with proper initial conditions in laboratory setups of con- 
densed matter systems as well as high energy systems. 
The required setup could consist of a high energy par- 
ticle entering a region of space where the presence of a 
strong external field would properly modify the field po- 
tential responsible for the dynamics. Alternatively, an 
abrupt change of an external parameter of a system (eg 
pressure) may lead to abrupt modification of the effec- 
tive potential determining dynamics. Such experimental 
setups could reproduce the conditions considered in our 
numerical simulations leading to soliton formation. In 
addition this soliton formation mechanism may be rele- 
vant in the early universe leading to topological defect 
formation at temperatures below the phase transition in 
regions of space where strong fields or pressure gradients 
are present. 

It is straightforward to extend our results to higher di- 
mensional systems. For example oscillation modes in sys- 



7 



v=0.6, (jj=0.3: No Breather Decay into Kinks 



v=0.5, ai=0.15: Breather Decay into Kinks 
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FIG. 8: a: The field evolution of the boosted breather with (v,u>) = (0.6,0.3), in a potential that switches from sine-Gordon 
(left) to $ 4 (right) at xq — 0. The 'particle' evolves oscillating without decay to solitons for these parameter values, b: The 
field evolution of the boosted breather with (v,lu) = (0.5,0.15), in a potential that switches from sine-Gordon (left) to 5> 4 
(right) at xo = 0. The 'particle' evolves oscillating but decays to a kink-antikink pair for these parameter values. 
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tems accepting vortex solutions may also decay to vortex- 
antivortex pairs for proper range of parameters v — lu in 
the context of a similar mechanism as the one described 
here. In fact, the reverse process of a vortex-antivortex 
annihilation to an oscillon was observed in Ref. [l7j |. 
The investigation of this type of oscillation mode decay 
to higher dimensional solitons is an interesting subject 
for future investigation. 

The Mathematica file used for the produc- 
tion of the figures may be downloaded from 



|http: / /leandros.physics. uoi.gr/partkinks.zip . 
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FIG. 9: The dots indicate parameter values where the decay of 
the gaussian oscillator into a kink-antikink pair was observed 
in the simulations at xo — where the potential switches to 
the $ 4 form. Dots indicate kink-antikink formation, with the 
different colours indicating the different potentials at which 
the decay takes place. Blue dots indicate decay in the region 
of sine-Gordon potential, red dots indicate decay in the region 
of the $ 4 potential, and purple dots indicate decay in both 
regions. Thick dots indicate formation of more complex prod- 
ucts. The cases where at least a kink-antikink pair is formed 
are indicated by a thick dot of the colour corresponding to the 
potential where the decay occurs, whereas the cases that only 
new oscillons are formed are indicated by pink think dots. 



Acknowledgements 

We thank T. Vachaspati, D. Steer and S. Nesseris for 
a careful reading of the paper which lead to insightful 
comments. C.S.C. thanks the Department of Physics at 
the University of Ioannina for hospitality during the de- 
velopment of this project. This work was supported by 
the European Research and Training Network MRTPN- 
CT-2006 035863-1 (UniverseNet). 



8 



[1] "Solitons and Instantons", R. Rajaraman, North- 
Holland, Amsterdam (1987). [10 

[2] "Kinks and Domain Walls", T. Vachaspati, Cambridge 

Universit y Press (2006). [11 

[3] D. Tong, |arXiv:hep-tIi/0509216l 

[4] S. Dutta, D. A. Steer and T. Vachaspa ti, Phys. Rev. Lett. [12 
101, 121601 (2008) [arXiv:0803.0670l [hep-th]]. 

[5] Y.M. Bunkov, H. Godfin (Eds.), 'Topological Defects [13 
and the Non-Equilibrium Dynamics of Symmetry Break- [14 
ing Phase Transitions', Kluwer Academic Publ., Dor- [15 
drecht /Boston/London 2000. 

[6] H. Arodz , J. Dziarmaga, W.H. Zurek (Eds.), 'Patterns [16 
of Symmetry Breaking', Kluwer Academic Publ., Dor- 
drecht/Boston/London 2003. [17 

[7] K. J. M. Moriarty, E. Myers and C. Rebbi, Phys. [18 
Lett. B 207, 411 (1988); M. Donaire, T. W. B. Kib- 
ble and A. Rajantie, N ew J. Phys. 9, 148 (2007) [19 
|arXiv:cond-mat/0409172]; W. H. Zurek Phys. Rept. [20 
276, 177 (1996) [arXiv:cond-mat/9607l35j ; V. M. Ru- [21 
utu, V. B. Eltsov, M. Krusius, Yu. G. Makhlin, B. Pla- 
cais and G. E. Vol ovik, Phys. Rev. Lett. 80, 1465 (1998); [22 
T. W. B. Kibble, |arXiv:cond-mat/0111082| I. Chuang, 
B. Yurke, A. N. Pargellis and N. Turok, Phys. Rev. E [23 
47, 3343 (1993). [24 

[8] M. B. Hindmar sh and T. W. B. K ibble, Rept. Prog. Phys. 
58, 477 (1995) [arXivihep ph/941 1342] . 

[9] R. Durrer, M. Kunz and A. Melchiorri, Phys. Rept. 364, 



H. 



Brandenberger, 



1 (2002) [arXiv:astro-ph/0110348] . 

L. Perivolaropoulos, Nucl. Phys. Proc. Suppl. 148, 128 

(2005) [arXiv:astro-ph/0501590l 
J. Magueijo and R. 
|arXiv:astrpph/0002030| 

D. G. Levkov and S. M. Sibiryakov Phys. Rev. D 71, 
025001 (2005) [arXiv:hep-th/0410198j . 

M. P. Mattis, Phys. Rept. 214 159 (1992). 

N. S. Manton and H. Merabet, |a~rXiv:hep-th /9605038 
T. Romanczukiewicz J. Phys. A 39, 3479 (2006) 
[arXiv:hep-th/0501066| . 

"Aspects of Symmetry: Selected Erice Lectures" , S. Cole- 
man, Cambridge University Press (1985). 
M. Gleiser and J. Thorarinson, |arXiv:hep-th/0701294| 
M. Gleiser, Phys . Rev. D 49; 2978 (1994) 
[arXiv:hep-ph/9308279] . 

W. Rindler, J. Denur, Am J. Phys., 56, (9), 795 (1988). 
http: / / www . wolf ra mxom / 1 

I. L. Bogolyubsky and V. G. Makhankov, JETP Lett. 24, 
12 (1976). 

E. J. Copeland, M. Gleiser and H. R. Mull er, Phys. Rev. 
D 52, 1920 ( 1995) [arXiv:hep-ph/9503217| . 

E. P. Honda, [a?Xiv:hep- ph/0009104p~ 

M. Hin dmarsh and P. Salmi Phys. Rev. D 74, 105005 

(2006) [arXiv:hep-th/0606016] . 



